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20,  Abstract  (Contd) 


bram  space  representation  allows  the  number  of  unknowns  to  be  reduced  from 
half  the  number  of  array  elements  to  the  number  of  Imposed  null  locations, 
and  results  in  a  sif^nificant  reduction  in  computation  time  from  that  required 
to  calculate  the  phase  perturbations  directly. 
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The  Use  of  a  Beam  Space  Representation  and 
Nonlinear  Programming  in  Phase^niy  Nulling 


1.  INTRODUCTION 

Thp  considerablo  rccont  interost  in  thf>  use  of  phasp-only  control  of  the  ele- 

1-9 

ment  weights  of  array  antennas  for  adaptive  nulling  reflects  the  growth  in 
importance  of  both  phased  arrays  and  adaptive  processing.  Phase-only  null 
synthesis  in  array  antenna  patterns  is  also  of  current  interest^®"*®  because  it 
can  help  establish  limits  to  what  can  be  achieved  adaptively,  and  possibly  suggest 
adaptive  procedures. 

Phase-only  null  synthesis  presents  analytic  and  computational  difficulties  not 
present  when  both  amplitude  and  phase  of  the  element  weights  can  be  freely  per¬ 
turbed.  The  restriction  of  the  weight  perturbations  to-be  of  the  phases  only, 
makes  the  nulling  problem  nonlinear  and  not  solvable  analytically,  Approximations 
and/or  numerical  techniques  must  be  used  to  calculate  the  phases  required  to 
impose  nulls  in  the  pattern.  As  with  nulling  with  combined  phase  and  amplitude 
control,  the  number  of  Imposed  nulls  is  typically  considerably  less  than  the  avail¬ 
able  number  of  degrees  of  freedom  (one  half  the  number  of  elements  for  phase-only 
nulling  in  real  antenna  patterns)  and  so  additional  conditions  must  be  Imposed  to 
define  a  unique  solution.  Since  in  either  null  synthesis  or  adaptive  nulling  it  is 


(necelved  for  publication  2  May  1983) 

Because  of  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 
See  References,  page  29. 
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generally  desired  to  preserve  the  original  pattern  Intact  as  much  as  possible 
apart  from  the  immediate  vicinity  of  the  imposed  null  locations,  a  natural  condi¬ 
tion  to  be  added  to  the  requirement  of  nulls  at  specified  location  '  is  that  the  per¬ 
turbation  of  the  pattern  be  minimized  in  a  mean  square  sense.  For  a  linear  array 
with  half-wavelength  intereleinent  spacing  this  is  equivalent  to  requiring  that  the 
sum  of  the  squared  absolute  values  of  the  weight  perturbations  be  a  minimum.  A 
more  general  condition  requiring  minimization  of  a  weighted  sum  of  the  squared 
absolute  weight  perturbations  can  provide  additional  flexibility  in  nulling  that  may 
be  desirable  in  certain  situations. 

In  a  series  of  recent  reports  we  have  investigated  various  approaches  to  cal¬ 
culating  the  phases  for  minimized  weight  perturbation,  phase-only  null  synthesis. 

A  straightforward  linear  approximation  method,  reasonably  effective  for  small 
phase  perturbations,  is  discussed  in  Reference  13;  a  method  based  on  iterated 
linearization  of  the  phase-only  nulling  equations  is  described  in  Reference  14;  a 
method  consisting  of  obtaining  the  best  phase -only  approximation  to  the  combined 
phase  and  amplitude  perturbation  solution  of  the  nulling  problem  is  reported  in 
Reference  17:  and  the  use  of  nonlinear  programming  techniques  to  calculate  the 
phase  perturbations  is  the  subject  of  Reference  18,  Of  these  methods  the  nonlinear 
progr:i  .ming  approach  is  the  most  general  and  effective.  In  this  report  we  de¬ 
scribe  a  variation  of  the  nonlinear  programming  method  based  on  the  work  of 
Baird  and  Rassweiler, 

In  their  basic  paper  on  phase-only  nulling,  ^  Baird  and  Rassweiler  considered 
the  problem  of  minimizing  the  mean  square  difference  between  a  desired  signal 
and  the  output  of  a  linear  array  by  varying  the  phases  only  of  the  array  elements. 
The  amplitudes  of  the  element  weights  were  assumed  equal  to  unity.  The  sources 
of  the  signal  and  the  interferences  were  modeled  as  discrete,  single -frequency 
directional  emitters.  The  vector  of  optimal  phases  was  shown  to  be  expressible 
as  the  phase  of  a  linear  combination  of  the  complex  conjugates  of  the  vectors 
giving  the  signal  and  Interferences  as  received  at  the  elements  of  the  array.  This 
representation  was  referred  to  as  the  "beam  space  decomposition,"  and  the 
coefficients  of  the  vectors  as  "beam  space  coefficients"  or  "beam  coefficients," 

The  purpose  of  this  report  is  to  apply  the  analytic  method  of  Baird  and 
Rassweiler  to  the  problem  of  Imposing  nulls  in  the  pattern  of  a  linear  array  of 
elements  by  varying  the  phases  only  of  the  element  weights,  subject  to  the  condi¬ 
tion  that  the  perturbations  of  the  weights  be  minimized  in  a  weighted  least-squares 
sense.  In  Section  2  we  obtain  a  representation  of  the  desired  phases  similar  to 
that  of  Baird  and  Rassweiler,  but  slightly  more  general  in  that  the  solution  allows 
for  an  arbitrary  amplitude  taper  of  the  element  weights,  and  for  any  choice  of  the 
coefficients  entering  into  the  weighted  sum  of  the  squared  weight  perturbations. 
When  the  phase  perturbations  are  small  the  perturbed  pattern  can  be  interpreted 
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as  the  original  pattern  plus  a  sum  of  cancellation  beams  directed  towards  the 
imposed  n-11  locations  »  id  the  locations  symmetric  with  respect  to  the  mainbeam. 
The  major  contribution  of  this  report  to  the  literature  on  phase-only  nulling  is 
contained  in  Section  3,  in  which  it  is  shown  that  a  highly  efficient  method  of  cal¬ 
culating  the  phase  perturbations  for  minimized  weight  perturbation  null  synthesis 
is  to  use  nonlinear  programming  computer  algorithms  to  calculate  the  coefficients 
in  the  beam  space  representation  of  the  phase  perturbations.  The  use  of  the  beam 
space  representation  enables  the  number  of  unknowns  to  be  reduced  from  half  the 
number  of  elements  to  the  number  of  imposed  null  locations.  Section  4  of  the 
report  contains  a  brief  discussion  of  the  use  of  the  beam  space  representation  in 
adaptive  phase-oniy  nulling. 


2.  ANALYSIS 

We  consider  a  linear  array  of  N  equispaced  isotropic  elements  (Figure  1) 
whose  field  pattern  is  given  by 

A  jd  u 
Po<“>  =  A  «n  ^ 

n- 1 


In  llq.  (1)  the  (a^j^are  the  complex  element  weights. 


N  -  1 


-  (n  -  1)  =  -d 


N-n+1 


n  =  1,  2, 


N 


(2) 


and 


u 


2y 


d  sin 


0  . 


where 

X  =  wavelength, 
d  =  interelement  spacing,  and 

0  =  angle  measured  from  broadside  to  the  array. 

The  phase  reference  is  taken  to  be  the  center  of  the  array.  Let  n  =  1,  2 . 

N,  be  the  set  of  perturbations  that  (a)  i-nposes  nulls  in  the  pattern  at  the  locations 
u  -  U|^,  k  =  1,  2 . K: 


N  N-l  N-2  3  2  1 


Figure  1.  Geometry  of  Array 


^  j*n 


(3) 


and  (b)  minimizes  the  weighted  sum  of  the  squares  of  the  absolute  values  of  the 
element  weight  perturbations 


n=l 


(4) 


The  weighting  coefficients  {c^}  in  Kq.  (4)  are  assumed  real  and  positive,* 

The  method  we  employ  to  derive  a  beam  space  representation  for  the  phase 
perturbations  is  that  of  Baird  a'  i  Rassweiler,  ^  We  let 


N 


(5) 


and  make  the  null  constraints  purely  real  by  multiplying  the  left-hand  side  of 
Eq.  (3)  by  its  complex  conjugate.  The  null  constraints  then  become 


N 


N 


EE 

m  =  l  n=l 


-j(d  -d  )u. 

m  n  m  n 


k  -  1,  2 . K 


(fi> 


The  choice  c^  "  1  for  all  n  makes  F  the  sum  of  the  squares  of  the  absolute 
values  of  the  weight  perturbations.  For  half-wavelength  spacing  of  the  arrav  ele¬ 
ments,  this  is  equivalent  to  minimizing  the  mean  square  pattern  pertui  bation  for 
9  from  -t/2  to  »/2.  Other  choices  of  the  c_  with  practical  application  to  null 

synthesis  are  also  possible:  for  example,  c  =  1/  a  2. 

n  '  '  n 


8 


•while 


=  E  cJaj2(2-2RelwJ)  . 


Since  the  {^jj}  acr  assumed  to  be  real  phase  perturbations,  the  {w^^}  by  Eq.  (5) 
have  unit  modulus.  This  is  expressed  by  an  additional  set  of  N  constraints 


w*Wn  =1  ,  n  =  1,  2,  N  . 


We  now  form  the  Lagrangian 


f'  • 


where  from  Eqs.  <fi)  and  (8)  respectively. 


N  N 

ES 

m  =  l  n-1 


a  a  e 
m  n 


w_  w_  ,  k  =  1,  2 . K 

m  n  ’ 


Dn  =  '*n'^n  *  ^  ^ . ^ 

The  {p|^)  and  {A^}  are  (real)  Lagrangian  multipliers.  A  necessary  condition  for 

F  to  have  a  stationary  point  is  that  the  gradient*  of  the  Lagrangian  with  respect  to 
T  19 

w  =  [Wj,  Wj . Wj^]  be  zero: 


^  20 
The  gradient  with  respect  to  a  complex  vector  z  is  defined  to  be 

^z  ^  ^Re(£j  ^^Im(zJ  * 

19.  Fletcher,  R.  (1981)  Practical  Methods  of  Optimization;  Vol.  2.  Constrained 

Optimization.  John  Wiley  &  Sons,  New  York,  Ch.  9. 

20.  Morse,  P.  M. ,  and  Feshback,  H.  (1953)  Methods  of  Theoretical  Physics. 

Part  I,  McGraw-Hill,  N.Y.,  p.  351. 
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(10) 


thrn  give* 


K 

2  5^! 

“  k»l 


•  2Rw 


(U) 


by  using  the  formula^ ^  for  the  complex  gradient  of  a  Hermltlan  quadratic  form. 
Finally, 


N 

\  '  2(A^Wj . 

~  n-1 


=  2  A  w 


where 


Substituting  f:qs.  (11).  (14).  and  (15).  in  Eq.  (10)  we  obtain 

i 

I 

.&  +  2Rw+2Aw*0  I 

or 

(R  +  A)  w  '  o  .  ^ 

i 

and  hence  j 


(R  +  A)~^a 


(15) 


(16) 


(17) 


Equations  (6).  (8).  and  (17)  form  a  system  of  2N  +  K  equations  for  the  2N  +  K 
unknowns  and  {p^|.  Since  the  constraints.  Eq.  (8).  specifying  the  use 

of  phase-only  weight  control  are  nonlinear,  the  system  of  equations  can  be  solved 

21.  Hudson.  J.  E.  (1981)  Adaptive  Array  Principles.  Institution  of  Electrical 
Engineers,  London,  p.  245. 
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only  by  th^  use  of  numrrlcal  techniques,  such  as  nonlinear  programming.  How¬ 
ever,  Eq.  (17)  can  be  used  to  derive  a  convenient  representation  of  the  phase 
perturbations.  Following  Baird  and  Rassweiler  we  define  a  sequence  of  matrices 

5m  “y 


5m 


m 


•kXkik 


».  2. 


k<l 


K  . 


5o  *  A  . 


with  the  vectors  defined  by  Eq.  (13).  Using  the  matrix  inversion  formula. 


22 


(A  +  a  z^)  =  A 


-1 


^ztA-‘z  ■ 


we  See  that  the  B  satisfy  the  backwards  recursion  relation 


«-l  .  T,-l  *'m5m-l-mJfm  2m^-l 

£m  •  £m-l  *“7T - TT-.T - 

V'  R  ‘v 


5o‘  ^  A‘‘ 


Since 


5k  *  A  +  R  . 


from  Eq.  (17) 


:  "  5k  £ 


(18) 


(19) 


B-X  ''K5K.liK2K5Kll 

-K-i  “7“ - rrri - 

^  *‘kiK5K-lXK 


22.  Householder,  A.S.  (1964)  The  Theory  of  Matric-s  In  Numerical  Anaivsis 
Blaisdell,  N.Y.,  p.  123: - - - ’ 
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£k-i  (- *  • 


where  Is  the  scalar 


,  vi  ,  o 


'K^K-K-l- 


We  now  repeat  the  process,  using  Rq.  (18)  and  letting 


Ik-i  ’2  •  ‘‘k^k  • 


so  that 


2  *  ?K-l^-l 


'‘K-l?K-2  5iK-liK.l5K-2 


’  £K-2{iK-l  *  ■>K.l^-l) 


2k-2|2  •  ‘•'K^K  "  '*^K-l^-l) 


where 


**K-l^-l-K-2  -iK-l 


Continuing  the  proc'-  ,  and  using  Rq.  (19)  at  the  end,  we  obtain 


2  -  A'M"  •  2  "^kik 

\  k»l 


fc3lwi*w<4v 


where 
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Assuming  that  the  Lagranglan  multipliers  of  the  constraint  functions  (0^1  of 
Kq.  (9)  are  all  positive,  the  factor  in  Eqs.  (21)  and  (22)  has  no  effect  on 

the  phases.  Hence,  the  phase  perturbations  are  given  by 


*n  '  ^n'^nl  “  **  "2 


.  n  =  I,  2 . N  .  (23) 


and  the  total  phases  of  the  element  weights  by 


+  t(  _ 

n  n 


-  phase  [c^a^  -  ^ 


,  n  •  1,  2 . N  .  (24) 


For  the  special  case  that  ja^l  "  Eq.  (24)  reduces  to  the  form  of  Eq.  (9) 

of  Baird  and  Hassweiler  (see  Appendix  B). 

The  coefficients  (>1^)  in  Eqs.  (23)  and  (24)  are  in  general  complex.  However 
in  the  important  special  case  that  (a)  the  unperturbed  pattern  is  real,  and  hence 
that  the  unperturbed  element  weights  are  conjugate  symmetric, 

®N-n+l  =  ,  n  =  1,  2.  ...,  N  ; 

and  (b)  the  weighting  coefficients  are  chosen  to  be  even  symmetric, 

‘^N-n+1  °  ®n  •  * . N  } 

23 

it  has  been  shown  that  the  phase  perturbations  are  odd -symmetric: 

^N-n+1  ””^n  '  n*l,  2,  N  .  (25) 

Since  the  {d^^}  In  Eq.  (23)  are  odd-oymmetrlc  by  Eq.  (2),  fo»  Eq.  (25)  to  hold  it 
is  then  necessary  that  the  coefficients  {yj^}  be  real. 

Equation  (23)  does  not  provide  an  explicit  solution  for  the  phase  perturbations 
since  the  coefficients  are  defined  in  terms  of  the  unknow^  Lagrangian  multi¬ 
pliers  Nevertheless,  the  form  of  the  phase  perturbations  la  usehil  because 

(a)  it  makes  possible  an  interpretation  of  phase-only  nulling  in|  terms  of  cancella¬ 
tion  beams;  (b)  it  can  serve  as  the  basis  of  a  numerical  method  for  calculating  the 


23,  Shore,  R.As  (1983)  A_  Pi^of  of  the  Odd^Symmetry  of  the  Pnases  for  Minimum 
Weight  Perturbation.  Phase-Only  Null  Synthesis,  RADC4TR-83~96. 
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phase  perturbations  for  null  synthesis  with  the  number  of  unknowns  equal  to  the 

number  of  imposed  nulls,  K,  rather  than  the  number  of  elements  In  the  array,  N; 

and  (c>  it  can  be  incorporated  into  adaptive  algorithms  for  phase-only  nulling. 

The  application  of  Eq.  (23)  to  null  synthesis  is  the  central  concern  of  this  report 

and  is  treated  in  Section  3.  The  use  of  Eq.  (23)  in  adaptive  nulling  is  briefly 

discussed  in  Section  4.  A  beam  space  interpretslion  of  phase-only  nulling  has 

been  given  by  Baird  and  Rassweiler,  and  also  in  Reference  13.  For  completeness, 

however,  we  give  it  briefly  here  as  well,  in  a  slightly  more  general  form. 

Assume  that  the  coefficients  are  real  and  that  the  unperturbed  phases, 

lij/  ).  are  of  the  form 
n ' 

'^n  '  ‘'^n“s  •  "  =  1.  2 . N  , 

so  that  the  main  beam  of  the  unperturbed  pattern  is  directed  towards  u  =  u^.  We 
rewrite  Eq.  (23)  as 


=  phase 


K 


‘=n'^J  -  2 

k  =  l 


=  tan 


r  K 

E 

>k  I'^n<“k  *  "s^ 

/  k=l 

K 

‘^nl^nl  • 

•E  ^k  ‘=o®td„(Uk-U3)l 

k=l 

n  =  1.  2 . N  .  (26) 


and  assume  that  the  phase  perturbations  are  small  so  that  can  l>c  approximateu 

by 


K 

S  ^k  l'ln<“k  -  “s>J 

k=l 


K 

=nl^ni  -2  T'k fV“k  ■  V' 

k  =  l 


,  n  —  1,  2,  N  , 


(27) 


and  the  weight  perturbations, 

(  e^*"  -  1  )  .  n  =  1,  2 . N  , 


n  n 
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approximated  by 


and  so  is  the  sum  of  K  pairs  of  beams,  onr  member  of  each  pair  directed  towards 
an  imposed  null  location,  u  =  Uj^,  and  the  other  member,  of  opposite  sign,  directed 
towards  the  symmetric  location  with  respect  to  the  direction  of  the  unperturbed 
main  beam. 
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u  =  -Uj^  +  2Ug  =  * 

The  shape  of  the  cancellation  beams  Is  determined  by  the  {cj^}.  For  c^  =  1, 

n  »  1,  2 . N,  and  the  magnitude  of  the  beam  coefficients  small  compared 

with  the  {la^jl  },  the  beams  are  of  the  form  sin  st"  ("^l  * 

corresponding  to  an  array  of  N  elements  with  uniform  amplitude.  For  other 
choices  of  the  {Cj^}  the  cancellation  beams  correspond  to  arrays  with  tapered 
amplitude  distributions.  Further  details  may  be  found  in  References  13  and  15. 

Although  the  interpretation  of  the  phase-only  nulling  cancellation  pattern  as 
a  sum  of  beam  pairs  is  based,  as  we  have  seen,  on  the  assumption  of  small  phase 
perturbations,  we  shall  for  convenience  in  this  report  refer  in  general  to  Eqs.  (23) 
and  (26)  as  "beam  space"  representations  of  phase-only  weight  perturbations,  and 
to  the  {7j^}  as  "beam  coefficients."  It  should  be  kept  in  mind,  however,  that  strictly 
speaking  the  beam  apace  interpretation  applies  only  to  phase  perturbations  small 
enough  for  the  approximations,  Eqs.  (27)  and  (28),  to  be  acceptable,  say 
4  I  <  0.  S  (rad).  It  should  also  be  noted  that  the  beam  space  representation, 

(29),  of  the  cancellation  pattern  is  in  general  a  nonlinear  superposition  of 
beams  because  the  beam  coefficients,  enter  into  the  expression,  Eq.  (30), 

for  the  Only  if  the  beam  coefficients  are  negligibly  small  compared  to  the 

{Cjjlajjl  }  does  the  representation,  (29),  become  a  linear  superposition  of 
beams. 

X  NUMERICAL  CALCULATION  OF  PHASE-ONLY  WEIGHT 
PERTURBATIONS  FOR  NULL  SYNTHESIS 

The  beam  space  representation,  (26),  obtained  in  the  previous  section  for 
the  phase  perturbations  for  minimized  weight  perturbation  null  synthesis,  can  be 
used  as  the  basis  of  a  highly  efficient  method  for  numerically  calculating  the  phase 
perturbations.  To  calculate  the  beam  coefficients,  (Tj^),  in  (26)  we  use  com¬ 
puter  algorithms  that  have  been  developed  for  solving  the  so-called  nonlinear  . 
programming  problem  —  the  problem  of  minimizing  or  maximizing  a  nonlinear 
ftinction  of  several  variables  subject  to  a  set  of  nonlinear  equality  and/or  inequality 
constraints.  Here,  the  nonlinear  hinction  we  wish  to  minimize  is  given  by  Eq.  (7), 
which  we  rewrite  in  the  form 
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(31) 


N 

n=l 


n  =  l 


and  the  nonlinear  equality  constraints  are,  from  Eq.  (3)  ai.d  the  odd-symmetry  of 
the  {djj}  and  the  {^^1.  ■ 

N 

=  S  l^nl  ‘=‘’®  I*n  ‘*n'“k  '  VI  "  ®  •  =  1-  2 . K  .  (32) 

n=l 

In  Reference  18  we  investigated  the  solution  of  this  non  inear  programming 

problem  when  the  U  )  are  the  unknowns.  The  performance  on  this  problem  of 

"  24  25 

the  two  nonlinear  programming  algorithms  LPNLP  and  VMCON  was  compared. 

Here  we  assume  the  form  of  (26)  for  the  phase  perturbations  with  the  K  beam 

coefficients  as  the  unknown  variables.  In  applications  to  null  synthesis  with 

large  arrays  and  the  number  of  imposed  nulls,  K,  small  relative  to  the  number  of 

independent  phases,  N/2,  this  approach  has  the  obvious  merit  of  reducing  the 

number  of  degrees  of  freedom  from  N/2  to  K.  Since  many  nonlinear  programming 

algorithms,  including  LPNLP  and  VMCON,  require  explicit  expressions  for  the 

derivatives  of  the  objective  function  and  of  the  constraint  Ainctions  with  respect  to 

the  unknowns  (some  algorithms  compute  these  derivatives  via  discrete  difference 

approximations),  we  give  expressions  for  the  derivatives  of  F  and  the  {EIjj}  with 

respect  to  the  in  Appendix  A. 

To  investigate  the  relative  performance  of  the  two  approaches,  oeam  space 
and  "phase  space,"  to  calculate  the  phase-only  weight  perturbations,  the  beam 
coefficient  method  was  used  to  calculate  the  phases  for  the  same  null  synthesis 
problems  as  were  used  in  the  study  described  in  Reference  18.  The  first  problem 
was  that  of  imposing  nulls  in  the  pattern  of  a  41-element  array  with  uniform  am¬ 
plitude  and  half-wavelength  interelement  spacing.  The  amplitudes  of  the  elements 


24.  Pierre,  D. A.,  and  Lowe,  M.J.  (1975)  Mathematical  Programming  Via 

Augmented  Lagrangians,  Addison -Wes ley,  Mass. 

25.  Crane,  R.  L.  et  al  (1930)  Solution  of  the  General  Nonlinear  Programming 

Problems  With  Subroutiy  VMCON.  Report  ANL-80-64.  ArgoMe  National 
Laboratory,  Argonne,  ill. 
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and  the  weighting  coefficients  were  set  equal  to  unity,  and  the  direction  of 

the  unperturbed  main  beam,  u^,  was  taken  to  be  zero  (that  is,  broadside  to  the 
array).  A  series  of  sets  of  imposed  null  locations  was  used  starting  with  one  null 
at  4.  0°,  then  two  nulls  at  4. 0°  and  4.  6°,  up  through  five  nulls  at  4. 0°.  4.  6°,  5.  2°, 
5.  8*^,  and  R,  4^.  Both  LPNLP  and  VMCON  were  run  in  double  precision  on  a 
CDC  RGOO  compute.'.  In  running  LPMLP  the  modified  Davidon -Fletcher -Powell 
(DFP)  conjugate  gradient  mode  (ISS  =  0)  was  used  without  automatic  reset  to  the 
gradient  direction  (IRESET  =  0).  The  convergence  parameters  were  set  at 

£,  =  1.0  X  lO”^®,  e,  -  f-i  -  1.  0  X  lO"*^.  In  running  VMCON,  the  tolerance  was 
^  -10  ^  ^ 

set  at  1.0  X  10  .  Both  programs  were  run  with  the  unknown  beam  coefficients 

set  initially  to  zero. 

In  Table  1  we  compare  null  depths  and  computation  times  obtained  for  this 
test  problem  with  LPNLP  and  VMCON  for  beam-space  and  phase-space  nulling. 
The  ph.ise-space  nulling  results  are  taken  directly  from  Reference  18.  The  most 
striking  result  is  the  remarkable  improvement  in  computation  time  obtained  with 
beam -space  nulling  as  compared  with  phase-space  nulling.  This  improvement  is 
especially  impressive  when  VMCON  was  used  as  the  nonlinear  programming 
algorithm:  for  the  one-,  two-,  and  three-null  cases,  computation  time  decreased 
by  a  fgetor  of  more  than  50,  while  for  the  four-null  case  the  computation  time  was 
reduced  by  a  factor  of  15.  Smaller  yet  considerable  decreases  in  computation 
time  were  obtained  when  LPNLP  was  used  as  the  algorithm.  Null  depths  obtained 
with  the  two  methods  are  comparable.  When  convergence  to  a  solution  occurred, 
the  beam  space  method  always  gave  the  same  phases  as  the  phase  space  method 
to  within  very  small  difference?.  The  five-null  case  could  not  be  solved  when  the 
beam  Space  method  was  used. 

To  help  understand  the  failure  of  the  beam  space  method  in  the  five-null  case, 
in  Tabie  2  we  have  tabulated  the  beam  coefficients  and  the  average  absolute  phase 
perturbation  as  a  function  of  the  number  of  imposed  nulls.  The  beam  coefficients 
are  tabulated  in  increasing  order  of  the  imposed  null  locations  towards  which  the 

i 

respective  cancellation  beams  are  directed.  The  average  of  the  magnitudes  of 
the  beam  coefficients  Increases  by  a  factor  of  more  than  70  and  the  average 
absolute  phase  perturbation  doubles  from  the  one-null  to  the  four-null  case,.  As 
discussed  in  Reference  18,  this  test  problem  was  intentionally  chosen  to  be  a 
difficult  one,  requiring  nulls  to  be  imposed  at  closely-spaced  locations  in  a  near-in 
sidelobe  region  of  a  uniform  array.  Viewed  intuitively  in  terms  of  a  picture  of 
cancellation  beams,  the  magnitudes  of  the  beam  coefficients  increase  strongly 
with  the  number  of  nulls  in  such  a  nulling  situation  because  the  main  lobes  of  the 
cancellation  beams  Interfere  with  each  other,  thus  requiring  extensive  mutual 
adjustment  of  the  magnitudes  of  the  beams  for  nulling  to  occur.  This  interference 
is  especially  evidenced  by  the  four-null  case  in  Table  2  in  which  not  only  are 
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Table  2.  Values  of  the  Beair  Coefficients  and  the  Average  Abso¬ 
lute  Phase  Perturbation  for  Beam  Space  Phase-Gnly  Nulling  in 
the  Pattern  of  a  41-Element  Uniform  Amplitude  Array  With  A/2 
Spacing.  Nulls  imposed  at  the  series  of  locations  4. 0  ,  4.  6  , 

5.  2°.  and  5.  8°  using  LPNLP  and  VMCON 


Number  of 
Nulls 

Beam 

Coefficients 

Average  Absolute 
Phase  Perturbation 
(rad) 

1 

-0.4 

0.28 

2 

-0.9 

0.31 

+0.5 

3 

+1.8 

-4.3 

0.33 

+2.  6 

4 

+  15.  6 
-42.9 
+42.0 
-14.9 

0.54 

the  magnitudes  of  the  beams  all  much  larger  than  that  required  to  impose  a  single 
null  at  the  respective  locations,  but  in  which  the  beam  coefficients  alternate  in 
sign  despite  the  fact  that  the  original  pattern  is  negative  at  the  locations  4. 0°,  4.  6°, 
and  5.2°,  and  is  positive  only  at  5.8°.  The  increase  in  the  magnitude  of  the  beam 
coefficients  with  the  number  of  closely-  ipaced  nulls  is  clearly  associated  with  an 
increase  in  numerical  difficulty  as  shown  by  the  increase  in  computation  time, 
especially  pronounced  in  going  from  the  three-  to  the  four-null  case.  The  numer¬ 
ical  difficulty  can  be  attributed  in  part  to  the  fact  that  Increasingly  large  changes 
in  the  beam  coefficients  are  required  to  bring  about  small  changes  in  the  objective 
and  constraint  functions.  However,  even  when  the  phases  are  taken  to  be  the 
unknowns,  the  five-null  cases  in  this  series  of  Imposed  nulls  presents  formidable 
difficulty  as  can  be  seen  in  Table  1  from  the  fact  that  only  LPNLP  converged 
to  a  solution  and  required  an  inordinately  large  amount  of  computation  time. 

Further  work  is  needed  to  develop  methods  for  handling  such  difficult  phase-only 
nulling  situations. 

The  second  test  problem  us«>d  to  compare  the  beam-space  and  phase-space 
nulling  methods  was  that  of  imposing  nulls  in  the  same  pattern  of  an  array  with 
41  elements,  half-wavelength  spacing,  and  uniform  amplitude,  at  the  series  of 


*lf  a  single  null  is  imposed  at  each  of  the  locations  4.0°.  4.6°,  5.2°,  and  5.8°, 
the  value  of  the  beam  coefficient  is  resp€-ctively  -0.41,  -0.31,  -0.14,  and  +0.07. 
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sets  of  locations  25^:  25°  and  35°:  and  so  on  up  through  a  set  of  five  imposed  nulls 
at  25°,  35°,  45°,  55°,  and  65°.  This  test  problem,  in  contrast  to  the  first  problem, 
was  chosen  to  be  relatively  easy,  requiring  nulls  to  be  imposed  at  wide-spaced 
locations  in  the  lower  sidelobe  region  of  the  pattern.  In  Tables  3  and  4  we  have 
tabulated  the  same  quantities  for  this  test  problem  as  in  Tables  1  and  2,  respec¬ 
tively.  Again  the  phase  space  results  are  taken  directly  from  Reference  18.  As 
was  noted  for  Table  1,  the  salient  feature  of  Table  3  is  the  remarkable  reduction 
in  computation  time  made  possible  by  employing  the  beam-space  nulling  method, 
especially  when  VMCON  is  used  as  the  nonlinear  programming  algorithm.  Even 
the  five-null  case  was  solved  by  VMCON  using  the  beam  space  method  in  about 
half  the  shortest  time  required  to  solve  the  one-null  case  when  the  phase  space 
method  was  used.  Examining  the  beam  coefficients  in  Table  4,  il  is  seen  that 
there  is  relatively  little  change  in  their  values  as  nulls  are  added:  the  coefficient 
associated  with  the  25°  beam  varies  between  0.067  and  0.065,  the  beam  coefficient 
associated  with  the  35°  beam  varies  between  -0,039  and  -0.  040,  and  so  on.  Again 
viewed  intuitively  in  terms  of  a  picture  of  cancellation  beams,  in  this  test  exam¬ 
ple  the  beams  are  coupled  only  through  their  sidelobes  instead  of  through  their 
main  beams  as  in  the  first  test  problem.  Hence,  the  cancellation  beams  are 
much  more  Independent  of  each  other  and  the  beam  coefficients  for  the  multiple 
Imposed  null  cases  are  very  close  to  the  values  they  have  when  only  a  single  null 
is  Imposed  at  each  of  the  respective  locations  (+0.067,  -0.042,  +0,055,  +0.029, 
and  +0.046  for  a  single  imposed  null  at  25°,  35°,  45°,  55°,  and  65°,  respectively). 
Because  of  the  relative  Independence  of  the  cancellation  beams,  little  mutual 
adjustment  of  the  beams  is  required  for  nulling  and  the  numerical  difficulty  of  the 
problem  Increases  only  slightly  as  additional  nulls  are  Imposed, 

It  may  be  of  Interest  to  note  that  the  values  of  the  beam  coefficients  for  the 
second  test  problem  can  be  estimated  rather  closely  from  a  simplified  form  of 
the  approximation,  (29),  for  the  cancellation  pattern.  In  this  example  the 
weighting  coefficients  (c^)  and  the  amplitudes  of  the  element  weights  { la^^l  )  are 
equal  to  unity,  the  unperturbed  main  beam  direction,  u  ,  equals  zero,  and  N,  the 
number  of  array  elements,  is  41.  Assuming  that  the  )  are  small  relative 

to  one,  we  can  approximate  the  {c^}  in  Eq.  (30)  by  c^^  *  1.  At  the  imposed  null 

locations  u  =  Uj^,  k  =  I,  2 . K.  the  cancellation  pattern,  Ap(u),  can  then  be 

approximated  by 

Ap(u^)  «  -  Y  7^  .  k  =  1.  2 . K  . 

if  we  neglect  the  contributions  to  the  cancellation  pattern  at  u  =  Uj^  of  the  beams 
directed  towards  the  other  imposed  null  locations  and  towards  the  locations 

23 


Table  4.  Values  of  the  Beam  Coefficients  and  the  Average  Abso¬ 
lute  Phase  Perturbation  for  Beam  Space  Phase-Otly  Nulling  In 
the  Pattern  of  a  41-Element  Uniform  Amplitude  Array  With  a/2 
Spacing.  Nulls  Imposed  at  the  series  of  locations  25°,  35°,  45°, 
55®,  and  (i5°  using  I-PNI.P  and  VMCON 


Average  Absolute 

Number  of 

Beam 

Phase  Perturbation 

Nulls 

Coefficients 

(rad) 

1 

+0. 067 

0.044 

2 

■t-0. 064 
-0.  039 

0.044 

3 

+0.066 
-0.  040 
+0.054 

0.053 

4 

+0.065 
-0.  040 
+0.052 
+0.018 

0.054 

5 

+0.065 
•0. 039 
+0. 052 
+0. 020 
+0. 046 

0.063 

symmetric  with  respect  to  the  main  beam  of  the  Imposed  null  locations.  Since  by 
definition  the  cancellation  pattern  equals  the  negative  of  the  value  of  the  unper¬ 
turbed  pattern  at  the  Imposed  null  locations, 

:ip{u^)  =  -Pq(u^)  .  k  *  1.  2,  ....  K  .  (33) 

It  follows  that 

>,j«-^Po<«g)  .  k  =  1.  2 . K  .  (34) 

Since  the  value  of  the  unperturbed  pattern  at  the  locations  25®,  35®,  45®,  55®,  and 
65®  Is  l.  41,  -0, 88,  1. 12,  0.  63,  and  0.  98,  respectively,  we  obtain  from  the  ap¬ 
proximation,  (34),  as  estimates  of  the  beam  coefficients  0.  069,  -0.043,  0.054, 

0.  031,  and  0.  048.  These  values  are  quite  close  to  the  calculated  values  given  In 
Table  4.  More  accura’**  estlTiates  can  be  obtained  from  the  approximation, 

(29),  by  replacing  the  {c^)  by  c^  as  above  but  not  neglecting  the  contributions 
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of  any  of  the  cancellation  beams.  Equation  (33)  then  yields  a  system  of  K  simul¬ 
taneous  linear  equations  that  can  be  solved  for  the  beam  coefficients. 


4.  ADAmVE  PHASE-ONLY  NULUNG 

The  application  to  adaptive  phase-only  nulling  of  the  beam  space  representa¬ 
tion  of  the  phase  perturbations  has  been  discussed  in  detail  by  Baird  and 
Rassweiler,  ^  and  v*e  will  here  touch  only  briefly  on  this  topic. 

In  the  previous  section  we  showed  that  computation  efficiency  for  minimized 
weight  perturbation,  phase-only  null  synthesis  in  purely  real  patterns  could  be 
significantly  Increased  by  using  a  beam  space  representation  of  the  phase  per¬ 
turbations,  thereby  reducing  the  number  of  unknowns  from  N/2  (N  the  number  of 
elements)  to  K,  the  number  of  imposed  nulls.  The  beam  coefficients  were  deter¬ 
mined  using  nonlinear  programming  algorithms  supplied  with  values  of  the 
objective  function,  the  constraint  functions,  and  their  derivatives  with  respect  to 
the  beam  ccefflclents.  In  adaptive  nulling  with  realistic  arrays,  the  pattern  can  no 
longer  be  assumed  real  and,  consequently,  the  beam  coefficients  in  the  beam 
space  representation,  Kq.  (23),  must  be  regarded  as  complex  constants. 

Equation  (23)  can  be  used  to  reduce  the  number  of  unknowns  from  N  to  2K  pro- 
vld*^!  that  the  directions  {Uj^{  of  the  sources  of  Interference  are  known  with 
sufficient  accuracy. 

In  adaptive  nulling  it  is  also  often  the  case  that  the  only  information  available 
to  the  adaptive  algorithm  (other  than  that  of  the  directions  of  signal  and  Jammers) 
is  that  of  the  total  output  power  of  the  array.  The  beam  coefficients  must  then  be 
determined  by  random  search  or  gradient  methods  that  attempt  to  reduce  the 
total  output  power  to  a  minimum.  One  danger  of  such  algorithms  is  that  the  de¬ 
sired  signal  will  be  nulled  out  in  the  course  of  reducing  the  total  output  power. 
Constrained  gradient  methods  such  as  that  developed  by  Frosi  are  not  helpfbl 
in  phase-only  adaptive  nulling  since  they  require  adjustments  of  the  element 
amplitudes  as  well  as  the  phases.  One  distinct  advantage  of  using  the  beam  space 
repr«'Sentation  for  phase-only  adaptive  nulling,  as  compared  with  adapting  with 
the  phases  themselves  as  the  control  variables,  is  that  in  certain  situations  (inter¬ 
ference  sources  spaced  not  too  closely  in  a  low  sldelobe  region  of  an  antenna 
pattern)  the  beam  coefficients  are  small  in  magnitude  and  the  cancellation  pattern 
then  consists  of  beams  with  complex  coefficients  directed  at  the  interference 
sources  and  at  the  locations  symmetric  with  respect  to  the  main  beam.  These 
cancellation  beams  can  have  only  a  small  effect  on  the  main  beam  and  so  the 

26,  Frost,  O.  L.  (1971)  An  algorithm  for  linearly  constrained  adaptive  array 
processing,  Proc.  IEEE  60:661-675, 
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retponae  of  the  antenna  to  the  desired  signal  is  not  likely  to  be  significantly  de¬ 
graded  as  a  result  of  the  adaptive  acQustment  of  the  phases.  The  choice  of  the 
weighting  coefflc'ents  (c^}  in  Eq.  <23)  can  be  used  to  control  the  shape  of  the 

cancellation  beams.  If  c  >1  for  all  elements,  the  cancellation  beams  are  of  the 

n 

form  sin  j /sin  j  corresponding  to  a  uniform  amplitude  array.  Other 
choices  of  the  (c  I  can  be  used  to  obtain  cancellation  beams  with  broader  main 

I  n  I 

beams  and  lower  sldelobes.  These  may  be  desirable  in  certain  applications. 
Further  details  may  be  found  in  Reference  13.  For  closely-spaced  jammers  the 
beam  space  method  is  likely  to  present  considerable  difficulties  in  adaptive  null¬ 
ing  because  the  beam  coefficients  are  then  large,  and  instabilities  can  easily  arise 
from  the  kind  of  phenomena  noted  in  the  previous  section  in  regard  to  the  first  test 
problem. 


S.  CONCLUSIONS 

In  this  report  we  have  describf'd  a  new  computational  method  for  minimized 
weight  perturbation,  phase-only  null  synthesis  in  linear  array  antenna  patterns, 
consisting  of  the  use  of  nonlinear  programming  computer  algorithms  to  calculate 
the  coefficients  of  a  beam  space  representation  of  the  optimal  phase  perturbations. 
The  beam  space  representation  enables  the  number  of  unknowns  to  be  reduced  from 
half  the  number  of  array  elements  to  the  number  of  imposed  null  locations.  The 
derivatives  of  the  objective  and  null  constraint  functions  with  respect  to  the  beam 
coefficients,  utilized  by  the  nonlinear  programming  algorithms,  are  calculated 
from  analytic  expressions.  Computation  time  is  significantly  reduced  in  general 
by  this  method,  compared  with  that  requlr<>d  to  calculate  the  phase  perturbations 
directly  with  nonlinear  programming  techniques. 
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Appendix  A 

Calculation  of  th«  Derivative*  of  the  Objective  Function  and 
Constraint  Function*  With  Respect  to  Beam  Coefficients 


In  this  appendix  we  obtain  expressions  for  the  derivatives  with  respect  to  the 
beam  coefficients  (7|^}  of  the  objective  function  F  given  by  Eq.  (31)  and  the  con¬ 
straint  functions  {E^}  given  by  Eq.  (32). 

From  Eq.  (31) 


^nl^nl^  ‘'’"^n‘4^  •  P  =  ^  2 . K  . 


and  from  Eq.  (32) 


^"^n 

~  ®^n  ®T'p 
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where 


I 

I 


a 
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k=l 


and 


K 

*n  ^  l‘^n<'>k  “  «s’> 

k=l 


Substituting  Eqs.  (A4)  and  (A8>  in  Eq.  (A6)  we  obtain 


coa‘  ^  ) 

^  =  — ;; - 2-  j  3ln  Idj^iUp  -  U^)l  +  cos  (d^iUp  -  u^)l  j 


^n 


(A9) 


The  derivatives  of  the  objective  function  and  constraint  functions  are  then  given, 
respectively,  by  Eqs,  (Al)  and  (A2)  with  -s- —  obtained  from  Eq.  (A9). 


Appendix  B 


Comparison  of  the  Minimized  Weight  Perturbation  Objective 
Function  With  the  Objectiva  Function  Used  by  Baird 

and  Rasmeiler 


As  noted  in  the  introduction,  Baird  and  Rassweiler^'*'  treat  the  problem  of 
minimizing,  by  phase-only  control,  the  mean  square  error  between  a  desired  sig¬ 
nal  and  the  array  output,  given  discrete  sinusoidal  sources  of  signal  and  interfer¬ 
ences,  In  the  special  case  that  the  Interferences  are  infinite  in  power,  the  mean 
square  error  can  be  minimized  only  by  placing  pattern  nulls  in  the  direction  of  the 
interference.  The  mean  square  difference  between  desired  signal  and  array  out¬ 
put  is  then  simply  |  p(u  )  -  p„(u„)I^,  where  p„(u„)  and  p(u„)  are  the  values  of  the 
original  and  perturbed  pattern,  respectively,  at  u  =  u^,  the  direction  of  the  desired 
signal.  For  the  main  beam  of  the  original  pattern  to  be  directed  towards  u  =  u  , 

the  unperturbed  weights  are  of  the  form  a  =  Ja  |  e  and 


|p(u,)  -p^(Us)|2  =  £ 


h  id  u 
n  s 
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Bl.  Baird,  C. A.,  and  Rassweiler,  G.G.  (1976)  Adaptive  sidelobe  nulling  using 
digitally  controlled  phase-shifters,  IEEE  Trans.  Antennas  Propag. 
AP-24:638-649.  - 
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where  the  are  the  phase  perturbations.  This  objective  function  differs  in 
general  from  the  minimized  weight  perturbation  objective  function  we  employ 
(Eq.  (4)J. 


n=l 


N 

=  2  2^  -coa*^)  .  {B2) 

n=l 

even  when  c^^  =  (a^^l  =1,  n  =  1,  ....  N,  the  case  treated  by  Baird  and  Rassweiler. 
Only  if  the  additional  restriction  of  the  phase  perturbations  to  be  odd-symmetric 
is  made  do  the  two  objective  functions  become  equivalent,  since  then 


Even  though  the  limiting  form,  for  infinite  power  interferences,  of  the  objec¬ 
tive  function  employed  by  Baird  and  Rassweiler,  Eq.  (Bl),  differs  from  the 
minimized  weight  perturbation  objective  function,  Eq.  (B2),  the  two  objective 
functions  lead  to  the  same  form  of  beam  space  representation  for  the  phase  per¬ 
turbations  when  c^  =  |a^|  =  1,  as  we  have  noted  in  Section  2  regarding  Eq.  (24). 
To  understand  why  this  is  so,  we  briefly  trace  through  the  derivation  of  Eq.  (24), 
using  Eq.  (Bl)  as  the  objective  function.  Starting  with  Eq.  (10),  the  equation  re¬ 
quiring  the  Lagranglan  to  have  a  zero  gradient,  it  is  seen  that  if  the  objective 
function,  Eq.  (Bl),  is  used  instead  of  F,  all  terms  on  the  right-hand  side  of 
Eq.  (10)  remain  the  same  except  for  which  is  replaced  by  where 

l.„l  (.’*"->)] 
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=  2  2  •“ml  l«nl<*m  - 
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using  Eq.  (5).  To  calculate  V^G  l^i 

w'  =  w  •*  _!.  • 

where 

1  =  (1.  1.  ...  ll"^  . 

so  that 

G  =  w'^  J?  w'  . 

where  Pis  the  symmetric  matrix 
P  = 

with 

Uji . Usif  • 

Then 

=>  2  P  w' 

=  2  P  (w  -  U  . 

using  the  formula  for  the  complex  gradient  of  a  Hermitian  quadratic  form  as  in 
obtaining  Eq.  (14).  Equation  (17)  is  then  replaced  by 

w  =  (P+ R  +  A)'‘(P  1)  . 


The  recursive  matrix  inversion  procedure  used  to  obtain  EXj.  (20)  now  gives 


2  =  A‘^  . 

where  the  {vj^}  are  given  by  Eq.  (13)  and  the  {yj^}  are  scalar  coefficients. 

it  follows  that 


where 


f/. 

y'i  *  -:tr^ -  .  k  =  1,  2 . K  . 

tk'y'o 

Hence,  just  as  Ekj.  (23)  is  obtained  from  Eq.  (20),  the  phase  perturbations 
given  by 


•  ph,..  2 


n  *  I,  2,  ....  N  . 


and  the  total  phases  of  the  element  weights  by 
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Since 


are 


(B3) 


Note  that  the  amplitudes  of  the  weights  do  not  appear  explicitly  in  these  expres¬ 
sions  for  the  phase  perturbations  and  total  phases  as  they  do  in  Eqs.  (23)  and  (24). 
The  values  of  the  coefficients  {7"  j  do,  however,  depend  on  the  weight  amplitudes. 
Comparing  Eq.  (B3)  with  Eq.  (24),  equating  with  -t  is  seen  that  for 

c  =  a  =  1.  the  two  solutions  are  identical  in  form.  Thus,  the  minimized  sig- 
nal  perturbation  objective  function,  Eq.  (Bl),  and  the  minimized  weight  perturba¬ 
tion  objective  function,  Eq.  (B2),  lead  to  the  same  form  of  representation  for  the 
phase  perturbations  when  c^  =  |a^|  =  1. 


